Insulating state and the importance of the spin-orbit coupling in Ca 3 CoRh06 
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We have carried out a comparative theoretical study of the electronic structure of the novel one- 
dimensional CasCoRhOe and CagFeRhOe systems. The insulating antiferromagnetic state for the 
CasFeRhOe can be well explained by band structure calculations with the closed shell high-spin d° 
(Fe 3+ ) and low-spin t% g (Rh 3+ ) configurations. We found for the CasCoRhOe that the Co has a 
strong tendency to be d 7 (Co 2+ ) rather than d 6 (Co 3+ ), and that there is an orbital degeneracy 
in the local Co electronic structure. We argue that it is the spin-orbit coupling which will lift 
this degeneracy thereby enabling local spin density approximation + Hubbard U (LSDA+U) band 
structure calculations to generate the band gap. We predict that the orbital contribution to the 
magnetic moment in CasCoRhOe is substantial, i.e. significantly larger than 1 fis per formula unit. 
Moreover, we propose a model for the contrasting intra-chain magnetism in both materials. 

PACS numbers: 71.20.-b, 71.70.-d, 75.25.+Z, 75.30.Gw 
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I. INTRODUCTION 

The synthesis of the one-dimensional cobaltate 
Ca3Co20gii2 and the discovery of the exotic stair-step 
jumps in the magnetization at regular intervals of the ap- 
plied magnetic fields have triggered a flurry of research 
activities, both experimentally^^^^^^^^i 6 - 
and theoretically . 17 ' 18 i 19 i 20 i 21 i 22 i 23 i 24 This compound con- 
sists of [Co2C>6]oo chains running along the c axis of the 
hexagonal unit cell, with in each chain alternating CoOg 
octahedra and CoOe trigonal prismsi The intra-chain 
coupling is ferromagnetic (FM) with a Curie tempera- 
ture of Tq = 24 K and the inter-chain antiferromagnetic 
(AF) with a Neel temperature of T N = 10 K£ Owing 
to an inter-chain magnetic frustration associated with 
the triangular lattice, the magnetic ground-state is ei- 
ther a partially disordered AF state or a spin-freezing 
stated Substitution of the Co by other transition met- 
als results in quite dramatic changes of the proper- 
ties. The CasCoRhOe compound has its Tq shifts up 
to a high value of 90 K and T N to 35 K^^2L2^30 
but the type of its magnetic ground-state is similar to 
Ca3Co2Q6 ) 26 ' 27 while Ca 3 FeRh06 is an intra-chain AF 
with Ttv = 12 Ki 25 i 29 ! 31 In these compounds the Rh ions 
occupy the octahedral sites while the Co/Fe reside within 
the trigonal prismsi 25 ' 29 

There is a debate about the electronic 
structure of these one-dimensional cobal- 
tates in the literature, in particular for 
Ca Co O J^- 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 ' 13 - 14 ' 15 - 16 ' 17 - 18 - 19 - 20 - 21 - 22 - 23 

and CaaCoRhOe^^^^^^^^ Consensus has 
yet to arrive about the valence, spin and orbital state 
of the Co ions, without which one could not make a 
reliable modeling of the magnetic properties. Focusing 
on the CasCoRhOe system, it was originally thought 
that the valence state is Co 2+ /Rh 4+ based on mag- 
netization measurements and bond valence sum. 25 ! 26 
Subsequent neutron studies concluded, however, that 
the Co 3+ /Rh 3+ scenario is a better description ! 27 i 30 
Very recently, a photocmission investigation favors yet 
the original Co 2+ /Rh 4+ assignment^ Existing band 



structure calculations also provide a mixed message: a 
generalizcd-gradient-approximation (GGA) study sug- 
gests the Co 3+ /Rh 3+ stated while an LSDA calculation 
leaves this issue open^ Important is to note that both 
calculations predict CasCoRhOe to be a metal, in strong 
disagreement with the experiment ] 10 ' 16 

To resolve these disagreements, we carried out a the- 
oretical study of the electronic structure and magnetic 
properties of CasCoRhOe in which we took the insu- 
lating state of the materia l 10 ' 16 as a key finding. For 
this we applied the LSDA+U metho d 33 ' 34 in order to 
take into account more explicitly the correlated motion 
of the electrons typical for transition metal oxides. Us- 
ing CaaFcRhOe as a reference, we find that the Co 
ions have a lower valence than Fe, and that the Co 
are in an orbitally degenerate state while the Fe are 
not. The valence assignment coming from our calcula- 
tions is the Fe 3 + /Rh 3 + for Ca 3 FeRh0 6 and Co 2 + /Rh 4+ 
for CasCoRhOg. The insulating state for Ca3FcRh06 
can be readily reproduced, but for CasCoRhOe we in- 
fer that the spin-orbit coupling (SOC) must be included 
in the LSDA+U scheme. We make a testable predic- 
tion, namely that there is a large orbital contribution 
to the magnetic moment in CasCoRhOe, about 1.7 [Ib 
per formula unit. We also calculate the various exchange 
constants and propose a model to explain the weak intra- 
chain AF in CasFeRhOe and enhanced intra-chain FM 
in Ca 3 CoRh0 6 . 



II. RESULTS AND DISCUSSION 

Our calculations are performed by using the full- 
potential augmented plane waves plus local orbital 
method^ We took the crystal structure data determined 
by x-ray diffraction^ The muffin-tin sphere radii are 
chosen to be 2.7, 2.0 and 1.6 Bohr for Ca, Co/Fe/Rh, 
and O atoms, respectively. The cut-off energy of 16 Ryd 
is used for plane wave expansion of interstitial wave func- 
tions, and 5x5x5 k-mesh for integrations over the Bril- 
louin zone, both of which ensuring a sufficient numeri- 
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FIG. 1: Density of states (DOS) of Ca 3 FeRh0 6 (panel a) and 
CaaCoRhOe (panel b) in the ferromagnetic (FM) state from 
LSDA. The Fermi level is set at zero energy. Both panels 
show, from top to bottom, the total DOS, the Rh 4d, Fe/Co 
3d, and O 2p partial DOS. Solid (dashed) lines depict the 
spin-up (down) states. 



cal accuracy. The SOC is included by second-variational 
method with the scalar relativistic wave functions^ and 
actually l z and s z are good quantum numbers due to the 
special trigonal crystal field of the Co sites as seen be- 
low. The easy magnetization direction is along the c-axis 
chains. In view of those c-axis chains being well sep- 
arated in the hexagonal a6-plane and of these two com- 
pounds behaving like a quasi-one-dimcnsional system, we 
study in this work the type and origin of the intra-chain 
magnetic coupling, thereby assuming a FM inter-chain 
coupling. 

As a reference, we first calculate the electronic struc- 
ture of Ca3FeRh06 in the FM state using the LSDA. 
Fig. 1(a) shows the density of states (DOS). One can 



clearly see that the octahedral Rh 4d has a large t 2g - 
e g crystal-field (CF) splitting of about 3 eVM With the 
t 2g shell completely full and the e g completely empty, 
the Rh is formally 3+ (4d 6 ) and takes the low-spin (LS, 
5=0) state. One can also observe that the Fe 3d spin-up 
states are essentially completely full and the spin-down 
completely empty. The Fe is thus formally 3+ (3d 5 ) and 
high-spin (HS, 5=5/2). The Fe partial DOS also shows 
that the unoccupied spin-down states are split roughly 
into two groups with about 1 eV separation. This is 
caused by the presence of a trigonal CF, making 
d-i) orbitals to lie higher than the nearly degenerate 
(d , d 2 , d_ 2 ) ,2ilZiiSil2>22i2ii22i2a Covalency reduces the Fe 
moment to 3.73 /ib but creates at the same time 0.20 
\xb moment at the Rh and 0.13 [Lb at each O, so that 
the total moment is still 5.00 fis, see Table I. The ma- 
terial is an insulator, since it is effectively a closed shell 
system due to the sufficiently large CF splitting at the 
Rh and exchange splitting on the Fe. We have also cal- 
culated the CaaFeRhOe in the intra-chain AF state, and 
found very similar DOS'ses (the spin part not consid- 
ered). Important is that the total energy of the AF is 
lower than the FM by about 11 meV (see Table I), ex- 
plaining why CaaFeRhOe is an AF. Our LSDA finding is 
in close agreement with the GGA results of Villcsuzanne 
and Whangbo.— 

Fig. 1(b) depicts the LSDA results for the Ca 3 CoRh0 6 
in the FM state. The DOS shows quite a number of sim- 
ilarities with that of the CaaFcRhOg, including the large 
CF at the Rh and the characteristic trigonal CF split- 
ting at the Fe/Co. The only significant difference is the 
fact that the Co 3d lies lower in energy than the Fe 3d, 
which is reasonable since Co is less electropositive than 
Fe. The consequences are, however, quite dramatic: the 
Fermi level now straddles through the (do, d 2 , d_ 2 ) part 
of the Co 3d spin down band, meaning that CaaCoRhOo 
would be a metal. This LSDA result is in clear disagree- 
ment with experimental observations.—^ Also Whangbo 
et aliU found a metallic solution. We have investigated 
the intra-chain AF state for CasCoRhOe and found that 
it is also metallic. It has, however, a higher energy than 
the FM by about 106 meV, see Table I. This means that 
system tends to be FM, at least if one is willing to trust 
these LSDA results as far as the magnetism is concerned. 
It is obvious, nevertheless, that this half-metallic solution 
is an artifact of the LSDA. 

The LSDA of Ca 3 CoRh0 6 shows that the Fermi level 
is located in a narrow Co 3d band with no more than 
1 eV width. It is then also natural to expect that al- 
ready modest electron correlation effects at the Co sites 
will have a pronounced effect on the energetics of this 
band. We therefore carried out LSDA+U band structure 
calculations^ for both the CaaCoRhOg and CaaFeRhOe 
with the Coulomb energy U is set to 5, 4, and 3 eV, and 
the Hund's exchange parameter Jh to 0.9, 0.9, 0.5 eV for 
the Co 3d, Fe 3d, and Rh 4d shells, respectively. 

We will first discuss the results for the CaaFeRhOg, our 
reference system. Fig. 2(a) shows the LSDA+U DOS of 
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FIG. 2: Density of states (DOS) of Ca 3 FeRh0 6 (panel a) and 
CaaCoRhOe (panel b) in the ferromagnetic (FM) state from 
LSDA+U. Inset: close up of the down-spin Co do, di and d-i 
DOS. See Fig. 1 caption for other notations. 



CasFeRhOe in the FM state. One can observe that the 
energy separation between the fully occupied tig shell 
and the fully unoccupied e g has increased in comparison 
with the LSDA results [Fig. 1(a)], in accordance with the 
inclusion of U. This also happens with the separation be- 
tween Fe 3d spin-up and spin-down states, whereby the 
spin-up states are pushed down to about -7 cV. This 
deep-lying state reflects the high stability typical for 3d 5 
ions, in which the high-spin configuration of the half-filled 
shell gives an energy gain of 4 Jh- With the Fe spin-up 
states being at the bottom of the valence band, the cova- 
lency with the Rh 4d and O 2p bands is much reduced. 
This is then also reflected in the increase (with respect 
to LSDA) of the Fe moment to 3.95 /is and decrease of 
the Rh to 0.11 [Ib and of each O to 0.11 /is- The total 
magnetic moment should not change: it remains indeed 



at 5.00 /j,b, see table I. As far as the low energy scale 
physics is concerned, the band gap has increased from 
about 0.2 eV to an appreciable 0.9 eV. It does not directly 
scale with the C/'s, since it is determined by the occupied 
Rh tig and the unoccupied Fe 3d spin-down states. Note 
that our conclusion that CasFeRhOg is an insulator with 
HS-Fc 3+ /LS-Rh 3+ and is a Hciscnberg spin-chain system 
(see below), is U- independent. 

Calculations for the intra-chain AF state of 
CasFeRhOe give very similar DOS'ses (the spin 
part not considered), with practically the same total 
energy as the FM (the difference being about 2 mcV is 
within the error bar, see Table I). This and the above 
LSDA result indicate that the intra-chain AF exchange 
interaction is very weak. In this situation the effects 
not included in our calculations, such as the detailed 
type of inter-chain ordering (we assumed that spins in 
all neighboring chains are in phase, i.e., we considered 
FM inter-chain ordering) may start to play a role and 
may finally determine the type of long-range ordering 
in Ca 3 FeRh0 6 . Thus from the present LSDA+U 
calculations we can only conclude that the FM and 
AF intra-chain orderings in this system are practically 
degenerate. We present below (Section III) qualitative 
arguments that in reality most probably the intra-chain 
interaction in CaaFcRhOg is weak but AF. All this is 
not inconsistent with the rather small Tjv value of 12 K. 

The LSDA+U results for CasCoRhOg are shown in 
Fig. 2(b). Similar to the CasFeRhOe case, the inclusion 
of U increases the splitting between the Rh 4d tig and e g 
orbitals as well as Co 3d spin-up and spin-down bands. 
Also similar is the reduction of the covalency of the 3d 
spin-up as it is shifted to very low energies. Surprisingly, 
however, the insulating state is not formed. The Fermi 
level still straddles through the lower-energy part of the 
Co 3d spin-down band. The main influence of U here is 
only that the effective crystal splitting with the (di, d_i) 
orbitals is increased. To investigate why there is no gap 
opening, we plot in the inset of Fig. 2(b) a close up of 
the down-spin Co do, di and d_2 states in the vicinity 
of the Fermi level. We can now observe that the spin- 
down do is fully occupied, and that it is the di and d_2 
bands which are partially occupied. This means first of 
all, with the spin-up states fully occupied, that the Co ion 
is essentially in the d 7 or 2+ valence state, and not in the 
d 6 or 3+ as proposed recentl y 17 i 27 i 30 Secondly, the inset 
reveals that it is the degeneracy of the di and d_2 bands 
which makes the LSDA+U to be inoperative to open the 
gap. If the do were higher in energy than the d2 and 
d_2, or if the Co valence were 3+, then LSDA+U would 
certainly have produced a band gap. For completeness, 
we have also calculated the AF state with the LSDA+U, 
and found also a metallic state with a 13 mcV higher 
total energy, see Table I. 

We propose to include the SOC to lift the di,d-i 
degeneracy. The use of the SOC to lift the degener- 
acy in correlated insulators has surprisingly been done 
only in a few instancesi 37 i 38 In most cases, this has 



been omitted, even for materials for which it is known 
that the orbital contribution to the magnetic moment 
is substantial .-^^^ii 1 .^ We claim that the SOC issue 
is essential for this particular CasCoRhOg compound, 
motivated also by the report that the magnetocrystallinc 
anisotropy is significant.— Note that for the Co 2+ ions in 
CasCoRhOe, the do and dijd-i states are split off from 
the higher-lying d\jd-\ states by about 1 cV [Fig. 1(b)], 
i.e., much larger than the SOC energy scale of about 
70 meVM Therefore, the SOC Hamiltonian can be sim- 
plified into just Q z s z by neglecting the and Z_s + 
mixing terms, both of which cause a mixing between the 
do and di/d_i, or between the c?2 (d-2) and d\ (d_i). 
Thus, the l z and s z are good quantum numbers in this 
particular case. 

The results of the LSDA+U+SOC calculations for 
CasCoRhOg in the intra-chain FM state are shown in 
Fig. 3. We now can observe that a gap has been opened 
near the Fermi level, consistent with the experimental 
finding that the system is an insulator . 10 ! 16 The band gap 
is small and is given mainly by the weakly split Rh Ad 
bands. The system is a Mott-Hubbard insulator in which 
the U in the Rh 4d shell is the determining factor^ As 
far as the Co 3d states are concerned, the band gap is 
about 2 eV, i.e. so large that it no longer determines the 
band gap of the system. A clearer view is offered in the 
inset of Fig. 3 which depicts a close up of the down-spin 
Co do, di and d-2 states in the vicinity of the Fermi level. 
One can see that the spin-down di band is now fully occu- 
pied and the d-2 fully unoccupied, i.e. there is no partial 
occupation anymore like in the LSDA+U. Interesting is 
to note that the spin-down do band has also an unoccu- 
pied component slightly above the Fermi level, but not 
at the Fermi level. This can be attributed to the band 
formation with the unoccupied Rh Ad 3z 2 — r 2 state (i.e. 
a\ g in the trigonal CF), demonstrating the Rh 4+ valence 
state with the t^g configuration and the one-dimensional 
character along the c-axis chain. 

While the approach to include the SOC is sound in 
view of the sizable magnetocrystalline anisotropy^ 6 - a 
definitive justification needs to come from a determi- 
nation of the magnitude of the orbital contribution to 
the magnetic moment: our LSDA+U+SOC calculations 
({7=5 cV for Co) predict for the Co 2+ ions not only a 
spin moment of 2.72 [Ib but especially a very large or- 
bital moment of 1.69 /ib along the c-axis chain. This 
is a testable prediction, which is actually confirmed by a 
recent x-ray magnetic circular dichroism study4£ We ver- 
ified that both the calculated spin and orbital moments 
stay constant within 0.1 /is when the U for Co is varied 
in the range of 2.5-6 cV. This huge orbital moment and 
the SOC firmly fix the orientation of the total magnetic 
moment along the c-axis chain, and thus the magnetism 
of CasCoRhOe is highly Ising-like. The calculations also 
predict an appreciable 0.54 [Ib spin moment on the Rh 4+ 
Ad, and 0.09 /is on each oxygen, showing the effect of 
covalency. Note that the Rh 4+ ions carry no orbital mo- 
ment, because the one hole resides on the a\ g singlet split 
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FIG. 3: Density of states (DOS) of Ca 3 CoRh0 6 in the ferro- 
magnetic (FM) state from LSDA+U+SOC. Inset: close up of 
the down-spin Co do, d2 and d_2 DOS. See Fig. 1 caption for 
other notations. 



off from the ti g manifold by the trigonal CF. We also note 
that the Co 2+ /Rh 4+ valence state of CasCoRhOg, as we 
obtained in this work, remain unchanged by counting 
the electron occupation of the down-spin Co do and di 
orbitals (the up-spin Co d shell is fully occupied) and the 
hole occupation of the down-spin Rh do (i.e. a\ g ) orbital, 
once an insulating gap is opened at U >2.5 eV for Co and 
U >2 eV for Rh in our LSDA+U calculations including 
the spin-orbit coupling. All in all, our Co 2+ /Rh 4+ so- 
lution is thus consistent with the large total magnetic 
moment as measured by neutron diffractioni 26 i 27 i 30 

We also have carried out LSDA+U+SOC calculations 
for CasCoRhOg in the intra-chain AF state, and we have 
found that this state is 70 meV higher in energy than the 
FM state, see Table I. This allows us to estimate the ef- 
fective FM exchange parameter of a Co 2+ pair being 15 
meV, using the simple Heisenberg model with the HS- 
Co 2+ spin-3/2. These numbers are substantially larger 
than for the Ca3Co20g system: there the AF-FM dif- 
ference is 12 meV and the exchange parameter about 
1.5 meV^ Our calculations thus nicely explain why the 
intra-chain Tq of CasCoRhOg (90 K) is much higher 
than that of Ca3Co20g (24 K). For completeness, we 
have also carried out LSDA+U+SOC calculations for the 
CasFeRhOg. The results, however, are almost identical 
to those of the LSDA+U, since the SOC is not operative 
for the closed spin-up Fe 3+ 3d shell and the full Rh 3+ t\ 
state. The absence of an orbital moment accounts for its 
Heisenberg spin character and very weak magnetocrys- 
talline anisotropy^i 
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TABLE I: Calculated electronic states of CasFeRhOe and CaaCoRhOe (FM: ferromagnetic; AF: antiferromagnetic; -Linsulator; 
-M:metal); total energy (meV) per formula unit relative to the ground state; spin moments (/is) at each Fe/Co, Rh, and O 
ion, as well as in the interstitial region and total magnetic moment per formula unit. Only the Co has the orbital moment and 
is shown. The Ca spin moment is about 0.01 [ib and the Fe orbital moment 0.02 /j,b, and are thus omitted. 
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III. MODEL OF MAGNETISM 

In order to explain both the experimentally observed 
and theoretically (computationally) confirmed weak in- 
trachain AF (Tn=12 K) in CaaFeRhOg and relatively 
strong FM (T c =90 K) in Ca 3 CoRh0 6 , we propose 
the following model. Starting from our finding that 
Ca 3 FeRh0 6 has the HS Fe 3+ and LS Rh 3 + (non- 
magnetic t^g) valence state, a charge fluctuation will cre- 
ate the Fe 2+ /Rh 4+ pair rather than Fe 4+ /Rh 2+ , since 
it is easier to stabilize a Rh 4+ rather than an Fe 4+ . In 
Fig. 4(a), we sketch the normal superexchange mecha- 
nism Fc 3+ -(0 2 ")-Rh 3+ -(0 2 ")-Fe 3+ via the intermedi- 
ate non-magnetic Rh 3+ (in the figure the O 2- is omitted 
for clarity). Since the HS Fe 3+ has a closed up-spin 3d 
shell, only the down-spin Rh 3+ electron [e.g., the a\ g 
(3z 2 -r 2 )] can make a virtual excitation, which naturally 
explains the intra-chain AF, with a low Tn due to the 
large Fe-Fe distance of 5.39 A . In fact, this result is a di- 
rect analog of the first Goodenough-Kanamori- Anderson 
rule4£ 




FIG. 4: Model of the intra-chain AF in Ca 3 FeRh0 6 (a) and 
of the intra-chain FM in CaaCoRhOe (b). 

CaaCoRhOg in contrast, has the HS Co 2+ and LS 
Rh 4+ (t|„, 5=1/2) valence state. We therefore con- 
sider the Co 3+ /Rh 3+ charge fluctuation (rather than the 
Co 1+ /Rh 5+ with a rarely observed Co 1+ state) for the 
virtual excitation path Co 2+ -(0 2_ )-Rh 4+ . Since Rh 4+ 
has a large t2 g -e g CF splitting which makes the empty 



e g states to be costly to reach, it will be the t\ g states 
which will allow for an electron transfer to the empty 
down-spin ai g orbital, as shown in Fig. 4(b). If the 
Co 2+ and Rh 4+ have a FM alignment, the down-spin 
Co 2+ dig (do) electron will do so; otherwise (in the AF 
Co 2+ /Rh 4+ alignment) the spin-up a\ g electron will do, 
but obviously with a higher energy cost than the former 
by 3 times Hund's exchange energy of the Co ion. Thus 
the Co 2+ and Rh 4+ strongly prefer the FM alignment, 
being an analog of the second Goodenough-Kanamori- 
Anderson rule4£ In addition to the spin-triplet oxygen- 
holc mediated intra-chain FM ordering present in both 
CasCoRhOe and Ca3Co20gr 3 . this mechanism accounts 
for a higher Tq in the former than in the latter. Note that 
if CaaCoRhOg had the HS Co 3+ and LS non-magnetic 
Rh 3+ valence states, as in the case of Ca3Co206 with 
the HS trigonal Co 3+ and LS non-magnetic octahedral 
Co 3+ , its Tq would be lower than that of Ca3Co206 due 
to a longer distance between the magnetic HS Co 3+ ions 
in the former than latter. 



IV. CONCLUSION 

In summary, we have performed systematic LSDA and 
LSDA+U calculations with the inclusion of the spin- 
orbit coupling for the new quasi-one-dimensional spin- 
chain materials CasFeRhOg and CasCoRhOg. We con- 
clude that CaaCoRhOg is a Mott insulator with the 
high-spin Co 2+ and low-spin Rh 4+ configurations, and 
that the correlated insulator CasFeRhOg has the closed- 
shell high-spin Fe 3+ and low-spin Rh 3+ . We predict 
that CaaCoRhOe has a very large orbital moment at the 
Co 2+ site, which also explains the strong magnetocrys- 
talline anisotropy and the highly Ising-spin like behavior. 
The inclusion of the spin-orbit coupling in the present 
LSDA+U calculations is also crucial to obtain the in- 
sulating state. This highlights the importance of the 
spin-orbit coupling in (nearly) degenerate correlated sys- 
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terns. In contrast, CaaFeRhOe has spin-only moments 
and behaves like a Heisenberg spin-chain system. More- 
over, our calculations reproduce the relatively strong fer- 
romagnetic intra-chain coupling in CaaCoRhOe but the 
weak, presumably antiferromagnetic one in CaaFeRhOfj, 
and we propose a model to explain the contrasting mag- 



netism. 
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